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In a scattering experiment, the induced dipole moments of a magnetoelectric point scatterer in response to 
driving fields are given by its polarizability tensor a. Its linewidth will be dictated by the local density of 
optical states (LDOS) at the scatterer' s position. To retrieve the magnetoelectric cross coupling components 
of a for an archetypical magnetoelectric scatterer — a split ring resonator — we study the frequency dependent 
extinction cross section cr ex t as a function of distance to an interface. Rather than following a purely electric 
or purely magnetic LDOS, we find a dependence which reflects the interplay of both dipole moments in a 
'mixed' magnetoelectric LDOS. For a strongly magnetoelectric cross-coupled microwave scatterer, we compare 
analytical point dipole with finite element method calculations. 



I. INTRODUCTION 

While electricity and magnetism are inextricably connected 
in optics, the interaction of light with matter is generally con- 
sidered to be almost entirely mediated by the photon's elec- 
tric field and the electric polarizability of matter 1 . In the last 
decade, this paradigm has shifted with the emergence of the 
field of metamaterial^H. In this field of research, complete 
control over the flow of light is promised by transformation 
optics 5 , provided one can engineer arbitrary AC (optical fre- 
quency oj) permittivity e{uS) and permeability p{uS) of the 
medium it passes through. To reach this goal, many work- 
ers nanostructure materials that intrinsically have p = 1 and 
e ^ 1 to create effective media that spoof a magnetic re- 
sponse p. The workhorse object in the field is the so called 
split ring resonator 6 16 , a metallic ring of an overall size of 
around A/ 10 with a single cut that provides a magnetic re- 
sponse through a circulating charge mode that corresponds to 
an LC-resonance. At telecom frequencies, experiments indi- 
cate that split rings, and similar metamaterial building blocks, 
have a strong magnetic polarizability of several times their 
physical volume^UM^. Therefore, it is possible to induce a 
strong magnetic dipole in such scatterers upon driving with 
the incident magnetic field of light. At the same time, meta- 
material scatterers often feature a strong electric, and a so- 
called 'magnetoelectric' polarizability, whereby electric driv- 
ing sets up a strong magnetic response and vice versa^. 

Parallel to the development of metamaterials, interest 
has recently emerged in engineering magnetic fluorescent 
transitions^" 2 ^. Just as the response of materials to light is 
dominated in nature by e, the fluorescence of deeply sub- 
wavelength objects like atoms, quantum dots and molecules, 
is usually entirely dominated by electric dipole transitions^. 
Thus, researchers in the field of spontaneous emission con- 
trol conventionally take the local density of optical states 
(LDOS), that quantifies how many photon states are avail- 
able for an emitter to decay into, as strictly meaning the lo- 
cal density of electric field vacuum fluctuations. This elec- 
tric field LDOS not only governs spontaneous emission, but 
is also commonly associ ated to, e.g., light generation by 
cathodoluminescencJSES, or the radiative damping of plas- 
monic, i.e., purely electrically polarizable, scatterers as mea- 
sured by Buchler et alPS. However, recent work by Taminiau 



et alphas revealed that in some rare earth elements magnetic 
transition dipoles can be sufficiently strong to be cleanly ob- 
served. Those transitions sample a different LDOS, namely, 
the local density of magnetic field vacuum fluctuations. 

In this Paper, the two developments described above come 
together in a single question and its answer. If a scatterer such 
as a split ring is indeed a magnetic, or even a magnetoelectric 
scatterer of mixed electric-magnetic character, which LDOS 
actually sets the radiative linewidth? To formalize this ques- 
tion, we ask how radiation damping affects a split ring if we 
abstract it to a point scatterer with a formally 6x6 polariz- 
ability^ of the form 



Q - {U Z) (5) ■ 



Here, the electric response to electric fields and magnetic re- 
sponse to magnetic fields are given by the 3 x 3 tensors ol e and 
glh, respectively. The off-diagonal quantifies how strongly 
a magnetic dipole can be set up by an electric field and vice 
versa. A polarizability such as Eq. ([T]) contains very nontrivial 
features, such as optical activity, pseudochirality, and handed 
radiation patterns, depe nding on the amount of magnetoelec- 
tric cross-coupling a J 17 * 27 * 1 In how far this polarizability 
truly describes experiments is a matter of debate. Matching of 
far-field transmission spectra of periodic arrangements of such 
metamaterial scatterers to a lattice model for point dipoles is 
excellent 29 but one may wonder if the dipole picture stands 
up to scrutiny in a near-field experiment. A particular near- 
field experiment would be to test if a split ring responds to the 
LDOS, a quantity specific to dipole transitions and scattering. 
In this paper, we first answer the question how a split ring's 
radiative linewidth is modified by the LDOS, and show that in 
addition to the electric and magnetic LDOS, a third quantity 
emerges in the form of a magnetoelectric LDOS. Secondly, 
we show that controlled variations in LDOS should allow one 
to measure the magnitude and test the conceptual validity of 
the point scatterer' s polarizability. Finally, we benchmark our 
proposal to finite element calculations. 
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FIG. 1 . Split ring resonators with two distinct orientations are placed 
above a perfect mirror. In the point dipole picture each split ring is 
described by an electric dipole p along the split ring gap (blue) and 
a magnetic dipole m pointing out of the split ring plane (red). The 
mirror images for both split ring orientations are depicted together 
with their respective image dipoles p' and m! . 
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II. MAGNETIC LDOS NEAR METALLIC AND 
DIELECTRIC INTERFACES 

The effect of the LDOS on a point scatterer is well un- 
derstood by first considering a polarizable dipole in front of 
a perfect mirror. As first demonstrated in a groundbreaking 
experiment by Buchler et al. 25 , the scattering resonance of a 
plasmon particle can be modified in width by mechanically 
approaching a planar reflective substrate. This effect can be 
interpreted in exactly the same manner as the explanation usu- 
ally given to Drexhage's experiment on the radiative transition 
rate of a fluorophore near a mirroi^Sl. The scattering reso- 
nance carries an electric dipole that hybridizes with its mirror 
image, which for a dipole moment parallel (perpendicular) to 
a perfect mirror interface has reverse (identical) orientation 
according to image charge analysis (Fig. [TJ. The two elec- 
tric dipoles together correspond to either a subradiant or su- 
perradiant configuration, depending on dipole orientation and 
distance. Consequently, the radiative linewidth oscillates with 
distance to the mirror in proportion to the electric LDOS. It is 
not immediately obvious that this method can be useful to also 
probe linewidth changes in objects with both an electric and 
magnetic dipole moment. While the rule for choosing the im- 
age dipole orientation reverses in the magnetic case, compared 
to the electric case, one should also consider that in a split 
ring the electric and magnetic dipole are at 90° relative orien- 
tation. In most experiments, the magnetic dipole is perpendic- 
ular to the substrate (labelled z-oriented from hereon), while 
the electric dipole is in-plane (x-oriented). An image dipole 
analysis assuming a perfect mirror predicts essentially no dis- 
cernible difference between the linewidth of an in-plane elec- 
tric dipole and an out-of-plane magnetic dipole. That electric 
and magnetic LDOS have essentially the same dependence for 
the perpendicular dipole orientations in a mirror charge anal- 
ysis was also noted by Karaveli et alP^. Discerning magnetic 
and electric LDOS effects hence either implies that one rotates 
the scatterer (Fig. [T] right) to have both electric and magnetic 
dipoles in plane, or requires that one finds an LDOS with dis- 
tinct electric and magnetic spatial dependence for the desired 
polarizations. 



FIG. 2. Relative electric (x-oriented dipole p, blue) and magnetic 
(z-oriented dipole m, red) LDOS for distance d from a vacuum/Ag 
(dashed lines) and vacuum/Si (solid lines) interface at a vacuum 
wavelength of 1.5 |^m. Especially at close distances d < 0.2 [im, a 
vacuum/Si interface provides a higher LDOS contrast of x-oriented 
electric dipoles and z-oriented magnetic dipoles than a vacuum/Ag 
interface. 



In Fig. |2j we plot the electric and magnetic LDOS for 
vacuum/Si and vacuum/Ag interfaces for both parallel and 
perpendicular dipole orientations with the aim of obtaining 
a large difference between the LDOS for x-oriented electric 
dipoles and z-oriented magnetic dipoles. To specify the cal- 
culation method for the LDOS at a position r above the in- 
terface, we calculated the imaginary part of the Green func- 
tion G(r, r), as described by Novotny and HechPl using the 
complex wave vector integration contour of Paulus et alP3. 
From hereon we suppress the argument of the Green func- 
tion. We generalize the calculation to encompass the electric 
LDOS IirGee, the magnetic LDOS IirGhh and the crossed 
Green dyadic ImG##, as specified in the appendix. We take 
e = 12.11 for silicon and e = —121.53 + 3.10z for silver, as 
tabulated for the resonant wavelength of 1.5 (arrP^Z], typical 
for 200 nm x 200 nm split rings, and plot LDOS normalized 
to the LDOS in vacuum (Fig. [2]). For the vacuum/Ag inter- 
face we observe that the magnetic z-oriented and electric x- 
oriented LDOS are quite similar in magnitude, except within 
50 nm of the interface. As anticipated from the perfect mirror 
intuition, a dielectric interface is advantageous in providing a 
higher contrast of magnetic electric LDOS contrast compared 
to a metal. Continuity conditions on E\\ and H\\ ensure that 
the electric and magnetic LDOS are highly distinct. The range 
over which a large distinction remains away from the inter- 
face extends well into the regime beyond the first oscillations 
in LDOS at 200 nm, as shown in Fig. [2] Therefore, scanning 
the separation distance between split ring and interface allows 
to independently vary the electric and magnetic LDOS over a 
substantial range. 
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III. RADIATIVE LINEWIDTH NEAR AN INTERFACE 

Now we proceed to examine the radiative linewidth of a 
point scatterer described by a magnetoelectric polarizability, 
as proposed by Sersic et alPl which is based on the static po- 
larizability introduced by Garcia-Garcia et alP generalized to 
include radiation damping. The induced dipole moments of a 
scatterer in vacuum are entirely set by its full 6x6 dynamic 



polarizability tensor that is of the form Eq. jlj>. In a 
quasistatic description of the scatterer, one starts from an LC 
circuit to obtain a static polarizability af^* mat consists of a 
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where 7]e,Vh,Vc are rea l parameters simply set by geometry. 
For an ideal, infinitely thin SRR all other elements of the 6x6 
polarizability are zero. To end up with an energy conserving 
scatterer, a radiation damping term must be added: 
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(3) 



where G is the 6x6 Green function. In vacuum the correc- 
tion amounts to the usual radiation damping term ilmG = 
2/3ik 3 I, where k = u/c, and I the identity matrix, that is 
standard in the field of plasmonics. In front of the interface, 
however, the Green function is modified according to Fig. [2] 
For the effectively 2x2 polarizability of the split ring, the rele- 
vant ImG tensor is also only 2x2, containing on the diagonal 
only the electric LDOS for x-oriented dipoles ImG ee,xx, and 
the magnetic LDOS for z-oriented dipoles IirGhh^zz^ while 
the off diagonal contains IulGeh,xz- 

We obtain the radiative linewidth as one would measure it in 
an extinction experiment, by calculating the extinction cross 
section in the following manner. The scatterer is driven by a 
plane wave impinging from above, plus its Fresnel reflection 
due to the interface. We take the incident beam as normal to 
the interface with the electric field polarized along the gap. 
We calculate the work done per unit cycle on the scatterer via 



W = ( ReE ■ Re^ 

dt 



ReH • Re 



dm \ 
~dt) 



(4) 



Plots of the work show an oscillatory dependence with dis- 
tance of the split ring to the interface, due to two effects. 
Firstly, the driving field forms a standing wave. Secondly, 
the polarizability varies with the oscillating LDOS. To obtain 
a true extinction cross section we divide out the local field 
strength of the driving 
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2Z 



(5) 



with Z the impedance of the host medium, in this case vac- 
uum. As Fig. [3ja, inset) shows, the retrieved extinction cross 



section is of Lorentzian spectral shape, and varies in width 
and central frequency as the scatterer is approached to the in- 
terface. We extract the resonance width, which is the sum of 
the radiative and absorptive damping rate of the scatterer. 

As a benchmark, Fig. [3] shows the linewidth of the extinc- 
tion cross section of a purely electric scattering sphere, i.e., 
taking t/e = 1, t\h = and r]c = 0, resonant at 1.5 |Ltm wave- 
length, with Ohmic damping rate 7 = 8.3- 10 13 s _1 and a par- 
ticle volume of V = 100 nm 3 . As the scatterer approaches 
the interface, its extinction linewidth oscillates, and almost 
doubles when the scatterer is close to the interface. We find 
that the linewidth T follows to great accuracy the dependence 
r = r a bs + r ra d x LDOS^^. Analogous to the calibration 
of quantum efficiencies of fluorophores^lES this dependence 
allows to extract the radiative and Ohmic damping rate of the 
particle, and thereby also the LDOS dependent albedo of the 
scatterer. For the sphere studied here, the albedo in absence of 
the interface is a = 0.41. This benchmark calculation shows 
that our model essentially reproduces the experimental obser- 
vation by Buchler et alH Similarly, a calculation for a purely 
magnetic scatterer verifies that the damping rate of a magnetic 
scatterer traces the magnetic LDOS [calculation not shown]. 

In Fig. [3j as a measure for T, we examine the extinction 
linewidth for objects that have both an electric and a magnetic 
character. For demonstration purposes, we take the electric 
and magnetic polarizability equally large at tje = Vh = 1- 
If no cross-coupling, i.e., no bianisotropy is present in the 
object (r]c = 0), the extinction linewidth simply traces the 
electric LDOS, provided excitation is normal to the sample 
so that the magnetic dipole is not directly driven at all [curve 
not shown]. As cross-coupling is introduced, and increased to 
its maximum value, the extinction linewidth shifts away from 
the purely electric LDOS, and towards the magnetic LDOS 
curves. At maximum cross coupling (rjc = 1), the extinc- 
tion linewidth exactly traces the mean of electric and mag- 
netic LDOS consistent with the fact that electric and magnetic 
dipole are equal in size. Generally, for this geometry and ex- 
citation, the averaging is weighted by the magnitudes of the 
dipole moments, i.e., \p\ 2 and |m| 2 . We conclude that an ex- 
periment such as performed by Buchler et al!^, but applied 
to a split ring can indeed provide a quantitative test of the 
magnetic and bianisotropic dipole response of a single object, 
and a calibration of the magnitude of polarizability tensor el- 
ements. 



IV. MAGNETOELECTRIC LDOS 

That a weighted average of electric and magnetic LDOS 
is obtained for a scatterer with both electric and magnetic 
dipole moment, may seem a likely general description of the 
physics describing radiative line widths. However, we note 
that in a split ring the electric and magnetic dipole are co- 
herently coupled through the magnetoelectric coupling, and 
hence may also probe the off-diagonal term in ImG. More- 
over, in a split ring the electric and magnetic response have 
a fixed phase relation, necessarily being a quarter cycle out 
of phase since the current inducing m and the charge sep- 
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FIG. 3. (a) Analytical calculation of the FWHM-linewidth of the extinction cross section normalized to driving field intensity as a function of 
distance d to a vacuum/Si interface for a purely electric (tje = 1, ?7iJ = rjc = 0) sphere resonant at 1.5 |^m wavelength, with Ohmic damping 
rate 7 = 8.3 • 10 13 s _1 and a particle volume of V = 100 nm 3 (dots) in comparison to the purely electric LDOS lineshape for this interface, 
(b) The dots represent the same quantitity as before, but for realistic split ring resonator (tjh — 0.7, t\e — 0.3, r\c — 0.4)^ that is oriented 
with the SRR plane parallel to the surface. Electric, magnetic and magnetoelectric LDOS are shown as lines, (c) Maximally cross-coupled 
split ring resonator (tje = tjh = T]c = 1) parallel to and with the gap pointing away and pointing towards the interface, respectively. 
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FIG. 4. FWHM-linewidth of the extinction cross section normalized 
to the local driving field intensity for a zero arm-length omega par- 
ticle oriented parallel to a Si interface as a function of distance d. 
Retrieved FWHM-line widths from FEM- simulation (dots) vs. point 
dipole calculation (lines). Insets show extinct power (left) as well as 
normalized extinction cross section (right) as a function of frequency 
/ and distance d as acquired from FEM-simulation. 



aration inducing p are related by charge conservation. This 
coherence for instance resul ts in a strongly handed response 
under certain viewing angle^^HD We predict that the coher- 
ence also affects the strength of interaction between the split 
ring and its mirror image in the substrate, i.e., the radiative 
linewidth. In the example we examined in Fig. [3] this ef- 
fect was fortuitously obscured due to the fact that the cross 
term Geh,xz by symmetry happens to be exactly zero. Mi- 
croscopically, the radiation emitted by an x-oriented electric 
dipole in front of a planar surface does not cause any mag- 
netic field along z to be reflected to the split ring to provide 
back action. If we rotate the split ring to stand with arms 
up, or point with arms down, cross coupling becomes impor- 



tant as Geh,xv 7^ 0- Microscopically, this indicates that an 
x-oriented electric dipole will receive a ^/-oriented magnetic 
field as parts of its reflection in the interface, which in turn 
will drive the object if it has a y-oriented magnetic polariz- 
ability. In Fig. [3] we report the radiative linewidth for a point 
scatterer oriented to exactly represent this case, i.e., that of 
a split ring that stands up with its legs normal to the plane. 
The radiative linewidth in this case does not trace a weighted 
combination of magnetic and electric LDOS. It will depend 
on Geh,xv instead. Remarkably, this dependence is differ- 
ent for the split ring pointing upwards to the split ring point- 
ing downwards although the object has the same electric and 
magnetic polarizability. The only difference is the sign of the 
cross coupling polarizability, i.e., whether the quarter cycle 
phase difference between magnetic and electric response is a 
lag or an advance. Thus the fact that radiation reaction is a co- 
herent effect means that the linewidth provides a direct way to 
measure phase relations between polarizabilities, and not only 
absolute values. For instance, one could measure if the quarter 
wave phase difference between cxe^cxh on the one hand and 
olc on the other hand, that is generally surmised from qua- 
sistatic ciruit theory for split rings, in fact carries over to real 
scatterers that are not negligibly small compared to the wave- 
length and that are not composed of ideal conductors. To our 
knowledge, this is the first proposition that a new property of 
the structure that can potentially be engineered independently 
of the well-known electric LDOS and the recently evidenced 
magnetic LDOS, may enter radiative linewidth modifications 
for dipole objects. We term this a magnetoelectric LDOS ef- 
fect. 



V. FINITE ELEMENT MODELLING EXAMPLE 

Our predictions for the effect of the magnetic and magne- 
toelectric LDOS on the linewidth of metamaterial scatterers 
are all subject to the assumption that the scatterers can ac- 
tually be described in a point dipole picture. An important 
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question that is yet to be tested in experiment and simulations 
is if this assumption holds at all for metamaterial scatterers, 
and if so for what classes of optical experiments. Therefore, 
we perform a numerical experiment and compare the radiative 
linewidth found from finite element method calculation with 
the point dipole predictions. To optimally discriminate for 
magnetic LDOS effects, we choose an omega-shaped parti- 
cle that earlier surface integral equation simulations(SIlP^) 
predict to have a very large magnetic and magnetoelectric re- 
sponse. In Ref!^, the polarizability was quantitatively re- 
trieved by projecting calculated scattered fields for the scat- 
terer held in free space on vector spherical harmonics. We 
have performed full-field finite element calculations of the 
omega- shaped scatterer, which is resonant in the microwave 
regime, again above a high-index substrate (n = 3.5). The 
scatterer geometry is that of a 60 nm thin flat loop of inner ra- 
dius 0.74 |Ltm and outer radius 1.19 |Ltm radius. Before the loop 
closes, the arms smoothly curve to be parallel over a length of 
390 nm, leaving a gap of 520 nm across. As material we use 
a Drude model for gold (e(uj) = — ujp/(uj(uj + ry)) with 
e b = 9.54, uj p = 2.148 • 10 15 s" 1 and 7 = 0.0092u; p . We 
employ the commercial COMSOL 3D FEM solver with ele- 
ments of quadratic order and a grid finesse down to 5 nm. Per- 
fectly matched layers enclose a cylindrical simulation domain 
(cylinder axis normal to the substrate) that extends 10 |Ltm 
around the scatterer. We perform a total field- scattered field 
simulation where, as in the point dipole model, the scatterer 
is excited by the superposition of a plane wave and its Fresnel 
reflection. We extract extinct power as the sum of scattered 
power (obtained from a near-field flux integral over a surface 
enclosing the particle) and absorbed power, and normalize ex- 
tinct power to the local driving strength. This procedure was 
tested on Mie scatterers to give cross sections to better than 
1%. 

Extinction spectra show a Lorentzian linewidth at all sep- 
arations, with a varying width and a slightly varying center 
frequency around a wavelength of 13.5 |Ltm. The center fre- 
quency varies because Eq. [3] in its most general form also 
contains ReG, corresponding to a real frequency shift due 
to the hybridization of the scatterer with its mirror image. 
We focus on the linewidth, plotted in Fig. [4] Evidently, the 
linewidth shows oscillations increasing in amplitude when ap- 
proaching the interface, and is twofold larger close to the in- 
terface than away from it. To move beyond this qualitative 
resemblance with the point dipole prediction, we also plot 
the linewidth found from point dipole calculations. No ad- 
justable parameters are used for the comparison, as we insert 
the polarizability values extracted from SIE calculations in 
Ref. [28l which are characterized through an/o^E = 0.3511 
and olcIole = 0.596. We note that the object has an on- 
resonance electric polarizability \ole\ = 3.9 |Ltm 3 , approxi- 
mately 30 times the particle volume. The point dipole model 
is seen to satisfactorily agree with the linewidth in simula- 
tions. Thereby, we establish that the point dipole approach not 
only describes far-field measurements on arrays of split rings, 
but is also directly applicable for split rings in the near-field 
of structures that modify LDOS. We propose that the resid- 
ual deviations contain interesting physics to pursue. Firstly, 



whether a metamaterial scatterer actually traces the magnetic 
and magnetoelectric LDOS can be seen as a fundamental test 
in the discussion in how far a spoof magnetic scatterer is ac- 
tually a true magnetic scatterer. Secondly, if one accepts that 
a scatterer that largely radiates according to the magnetic and 
magnetoelectric LDOS is a bona fide magnetic dipole, one 
can assess on basis of the residual deviations between simu- 
lation and analytical model in how far multipolar corrections 
are important. 



VI. CONCLUSION 

To conclude we have examined the dependence of the ra- 
diative linewith of split ring scatterers on their distance to an 
interface that modifies the electric LDOS, the magnetic LDOS 
and a new quantity that we term magnetoelectric LDOS. We 
propose that this linewidth, i.e., the backaction of the field 
radiated by the scatterer on itself, can serve as a calibration 
probe of the complex polarizability tensor and as a funda- 
mental test of the proposed dipolar nature of metamaterial 
scatterers. Of particular note is the concept of magnetoelec- 
tric LDOS, whereby the electric dipole of an object radiates 
magnetic field that back-acts on the magnetic dipole. It is an 
interesting question to explore whether such a magnetoelec- 
tric LDOS will also affect quantum mechanical transitions. 
While in the recent breakthroughs by Taminiau et al. and Kar- 
aveli e t al. ma gnetic-only transitions in rare earth ions are en- 
hancecPDSEl, it is an open question if transitions with a clear 
simultaneous electric and magnetic character can be found 
Conversely, we propose that coupling single emitters with a 
purely electric response to magnetoelectric scatterers may al- 
low to spoof quantum mechanical transitions with a magneto- 
electric character, as antennas tend to imbue their polarization 
characteristics on emitters. Such emitters woud likely have 
interesting chiral properties, since magnetoelectric cross cou- 
pling implies optical activit}^^. 
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Appendix: Green function 

To calculate LDOS, we require the 6 x 6 dyadic Green func- 
tion near a planar interface where source r' and observation 
point r are in the same half space. We separate the Green func- 
tion in a free part (in absence of an interface) and a reflected 
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part 



with 



G(r, r') = G free (r, r') + G re fl ec ted(r, r') 



G free (r,r)=^ ^ + yy j G(r, r ) (A.l) 

where k = uju/c is the wave number in the medium of index n 
that contains both r and r', c is the speed of light, and G(r, r') 
is the scalar Green funtion. The reflected part of the green 
function reads 

ik 2 f°° 

G reflected (r,r) = — J k\ | dk\ i [j (k\ i R)M + 

j 1 {k n R)M 1 +j 2 {k n R)M 2 ]e ik ' z 

(A.2) 



where if r = (x, z) and r' = (V, y' z') we define cylindri- 
cal coordinates through (i? cos 0, sin 0, Z) = (x — x f ,y — 
y' ,z + z'). With fcs we denote y//c 2 — ll^|| l| 2 > while j n (x) is 
the Bessel function of order n. The 6x6 matrices contain 
the fc|| dependent Fresnel reflection coefficients r s and r p for 
s and p polarization. In detail: 
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r s /k z -r p k z 
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and finally also 



Mo 



({rs/kz + r p /c 2 ) cos 20 (r s /fc^ + r p /c 2 ) sin 
(rs/k z + sin 20 


— (^a + ^p) sin 20 
(r s + ^p) cos 20 

v o 



(r s +r p ) sin 20 — (r s + r p ) cos 20 0\ 

(r s /k z + r p k z ) cos 2(j) — (r s + r p ) cos 20 — (r s + r p ) sin 20 



(r s +r p ) cos 20 (r s k z + r p /k z ) cos 20 (r s /c 2 + r p /k z ) sin 20 

(r s +r p ) sin 20 (r s /^ + r p /k z ) sin 20 -(r s k z + r p /k z ) cos 20 

0/ 



r 



Throughout we have used the units of RefP 
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